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Abstract —In this study, energy-efficient deterministic adap¬ 
tive beamforming algorithms are proposed for distributed sen¬ 
sor/relay networks. Specifically, DBSA, D-QESA, D-QESA-E, 
and a hybrid algorithm, hybrid-QESA, that combines the ben¬ 
efits of both deterministic and random adaptive beamforming 
algorithms, are proposed. Rigorous convergence analyses are 
provided for all our proposed algorithms and convergence to the 
global optimal solution is shown for all our proposed algorithms. 
Through extensive numerical simulations, we demonstrate that 
superior performance is achieved by our proposed DBSA and 
D-QESA over random adaptive beamforming algorithms for 
static channels. Surprisingly, D-QESA is also more robust against 
random node removal than random adaptive beamforming algo¬ 
rithms. Eor time-varying channels, hybrid-QESA indeed achieves 
the best performance since it combines the benefits of both types 
of adaptive beamforming algorithms. In summary, our proposed 
deterministic algorithms demonstrate superior performance both 
in terms of convergence time and robustness against channel and 
network uncertainties. 

Index Terms - Beamforming, convergence analysis, dis¬ 
tributed algorithms, time-varying channels 

I. Introduction 

There has been a demand for improving the reliability 
and energy efficiency of wireless communication systems 
pertaining to distributed sensor/relay networks. A lack of 
constant power supply in a distributed wireless system, such 
as a wireless sensor or a relay network, and the possibility 
of hardware failures have driven the emergence of energy- 
efficient algorithms for such applications. In particular, dis¬ 
tributed beamforming has been proposed as a viable solution 
where all distributed transmitters seek to align in phase at 
the receiver end. However, it is difficult to implement such 
transmit beamforming scheme in a distributed fashion in prac¬ 
tice because perfect channel state information (CSI) needs to 
be made available at all distributed transmitters. Tremendous 
overhead is required to feed back CSI from the receiver to all 
distributed transmitters. 
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In literature, the efforts of designing efficient distributed 
adaptive beamforming algorithms can be categorized into two 
classes: random adaptive beamforming algorithms or deter¬ 
ministic adaptive beamforming algorithms. The latter has been 
greatly investigated in the literature over the last decade. In ||T), 
a randomized approach was proposed in which the receiver 
uses a feedback link that allows each transmitter to make 
proper phase adjustment. In Q, Mudumbai et al investigated 
an adaptive distributed beamforming paradigm that requires 
only a single bit feedback. The scalability and the rate of 
convergence analysis of that scheme were given in By 
considering this algorithm as a local random search algorithm, 
Q provided a comprehensive analysis of the fast convergence 
of the algorithm. An extensive convergence analysis of these 
scheme was presented in Q. A new scheme inspired from 
a completely different field was proposed in ||^, where a 
bio-inspired robust adaptive random search algorithm was 
presented, and has been proven to converge in probability. 
Other examples of feedback-based synchronization procedures 
include Il7|-|[T§. 

Relative few efforts have been spent on the studies of 
deterministic adaptive beamforming algorithms. Several de¬ 
terministic algorithms for distributed beamforming have been 
proposed in the literature. For example, Thibault et al in¬ 
troduced a deterministic algorithm with individual power 
constraint 0- This algorithm was extended to the case of 
time-varying channels in (T^. Simulation results show that the 
deterministic approaches outperform the random methods in 
the case of static and time-varying channels. For amplify-and- 
forward wireless relay networks, an algorithm using additive 
deterministic perturbations was presented in 0. Fertl et al 
further investigated a multiplicative deterministic perturbations 
for distributed beamforming under a total power constraint 
|14|. Although faster convergence is often observed for de¬ 
terministic adaptive beamforming algorithms, they are more 
sensitive to channel uncertainties in general. 

In this study, we investigate the design and analysis of 
deterministic adaptive beamforming algorithms. We propose 
a Deterministic Bisection Search Algorithm (DBSA) that is 
inspired by the idea of bisection search. Furthermore, we 
propose a Deterministic Quadratic Equation Search Algorithm 
(D-QESA), inspired by a shift of perspective that views 
each transmission as a RSS function evaluation and views 
the problem of adaptive beamforming for each distributed 
transmitter as solving quadratic equations with independent 
variables. Eor the case of equal channel gains, we propose a 
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modified version of D-QESA, termed D-QESA-E, to speed up 
convergence for this setting. Note that deterministic algorithms 
usually have worse performance for time-varying channels. 
Therefore, we further propose a hybrid algorithm, hybrid- 
QESA, that combines the benehts of both deterministic and 
random adaptive beamforming algorithms. 

The remaining part of this paper is organized as follows. 
In Section II, we presents a system model for the adaptive 
distributed beamforming problem and formulate the phase 
alignment problem. The details of DBSA is described in Sec¬ 
tion III. In Section IV, we propose D-QESA, D-QESA-E, and 
hybrid-QESA. Convergence analysis is conducted for all our 
proposed algorithms in Section V. Extensive simulation results 
are presented in Section VI to demonstrate the performance 
advantages of our proposed algorithms. We conclude the paper 
in Section VII. 


we denote RSSmaa;[t2] as the global maximum value of the 
RSS function, i.e.. 


RSS, 




Ns 
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We then denote the beamforming gain ratio p[n] as 


Pn, [n] 


RSS(^i[n],-- - ,givjn]) 

RSS 

max [^] 


(4) 


The beamforming gain ratio represents how close to the 
global maximum the considered algorithm reaches, and attains 
its maximum value of 1 when perfect phase alignment is 
achieved. Note that RSSmaa:[?T-] and [n] can vary with time 
if we consider uncertainties in the wireless network topology 
or time-varying channel coefficients. 


II. System Model and Problem Eormulation 


In this study, we consider a wireless sensor/relay network 
consisting of distributed transmitters where a common mes¬ 
sage s S C is to be conveyed to the receiver end. Let Ng 
be the total number of distributed transmitters, each with an 
average power constraint of E'[|sp]. Each transmitter and the 
receiver is assumed to be equipped with single antenna. We 
assume that the channel between them is frequency flat and 
slow fading. In addition, there exists an error-free feedback 
link from the receiver to all distributed transmitters. 

According to these system specihcations, the discrete-time, 
complex baseband model is given by 

Ns 

uM = E ^iMgi[n]s + w[n] 

i=i 

= 2 + w[n\ 

i^l 

where y[n\ denotes the received signal, hi[n\ = 
corresponds to the channel coefficient, gi[n\ = 
corresponds to the beamforming coefficient, and w[n] is the 
additive white Gaussian noise for the f-th distributed transmit¬ 
ter. 

To simplify notation, let 9i [n] = [n] -I- tpi [n] be the total 

phase of received signal from f-th transmitter during the n-th 
transmission. Also, we assume that s = y/P and impose a 
hxed power constraint bi = 1 among transmitters. Therefore, 
the received signal can be expressed as 

Einally, we assume that the strength of the composite signal 
from distributed transmitters can be perfectly estimated at the 
destination and the received signal strength (RSS) function can 
be described as 


RSS(6»i[n],--- ,0Ns[n]) = '/P 


2=1 


( 2 ) 


Note that the primary goal here is to maximize the RSS 
function given in Q for distributed beamforming schemes to 
exploit the potential power gain efficiently. Eurthermore, it is 
clear that to reach the global maximum, all phases need to be 
completely aligned, i.e, 9i[n] = 02[n] = • • • = 9N^[n]. Here, 


III. Deterministic Biscetion Search Algorithm 

In this section, we introduce in details the proposed DBSA 
algorithm, which is inspired by the idea of bisection search. 
To simplify notation, let us assume that the phase of each dis¬ 
tributed transmitter remains conhned within an angle varying 
betweens [0, 27r]. By viewing the phase alignment problem as 
a maximization problem within an Ng dimensional hypercube 
with an edge length of 27r, it is then suggested that the 
searching process for the global maximum can be done in 
a bisection fashion. 

The phase alignment procedure is done in a greedy fashion 
within the proposed DBSA. In fact, only one distributed 
transmitter is permitted to adjust its phase at any time while 
all other transmitters keep their phase unchanged. With the 
pseudocode given in Algorithm 1, we elaborate the algorithm 
by providing a detailed step-by-step description as follows; 

Step 1- Initialization: The system is initialized 
by randomly generating the initial phase 0[O] = 

{0i[O],--- ,6>Arj0]} G The phase adjustment parameter 
a, known to all transmitters, should be configured to a 
number ainit that divides 27r. Distributed transmitters jointly 
beamform a hxed symbol s = y/P using the generated phase 
conhguration. The receiver then records the value of the RSS 
function given by corresponding to the initial phase 0[O]. 

Step 2- Initial Phase Rotation: The algorithm undergoes 
Ng rounds of phase evaluations corresponding to the Ng 
distributed transmitters. 

During the i-th round, only the i-th transmitter transmits 
for K — 1 iterations. The f-th transmitter alters its phase by 
iterating through the elements in the set 0^ given by: 

= {0*[O] + «,•••, 0,[O] + iK- l)a} Wi < Ng (5) 

where K is an integer that equals to ja. At each iteration, 
the receiver compares the RSS function with the largest 
RSS function recorded so far, and feedback a single bit that 
indicates whether the current selected phase conhguration by 
the i-th transmitter exceeds previous recorded RSS function. 
Should the RSS function exceed the largest recorded value, 
the feedback bit is set to 1, and to 0 otherwise. 

By the time all phase conhgurations in 0^. are tested, 
thei-th transmitter should be able to distinguish the phase 







3 


configuration that correspond to the best RSS performance, 
and this particular phase configuration should replace its initial 
phase configuration generated in Step 1, i.e., the i-th element 
of 0[O]. Once all transmitters are done with the evaluation 
procedures, the parameter a is halved, and the algorithm enters 
Step 3. 

Step 3- Forward Adjustment: Before proceeding to Step 
5, the algorithm undergoes iterations where each iteration 
involves procedures Step 3 and Step 4. At the n-th iteration, 
the n-th distributed transmitter adjusts its phase by adding 
up the parameter a to its phase, while all other distributed 
transmitters keep their phases unchanged. From the system 
point of view, the adjusted phase vector can be expressed as: 

e'[n] = [9i[n- I],--- ,6n-i[n- l],6»„[n- 1] + a, 
6i„+i[n- 1], - • • ,6»AfJn- 1]]'^ 

The receiver compares the RSS function corresponding to 
the adjusted phase vector 9 [n] with the highest recorded RSS 
function up to the current iteration, and broadcasts a one-bit 
feedback signal Bp back to all transmitters. If improvement 
on the RSS function is observed, the adjusted phase is kept 
by the n-th transmitter and the algorithm proceeds to the next 
iteration that adjusts the (n -|- l)-th transmitter provided that 
n < Ns- Otherwise, if no RSS improvement is observed, the 
algorithm enters Step 4, or if n = the algorithm enters 
Step 5. 

More precisely, the feedback signal Bp k generated as 

Bp = l{RSS(0'[n]) > max RSS(0[i])} n\ 

Step 4- Reverse Adjustment At this stage of the algorithm, 
the n-th distributed transmitter adjusts its phase by adding the 
inverse of the adjustment parameter to its phase, where the 
inversely adjusted phase vector 9"\n] can be expressed as 

e"[n] = [6»i[n - 1], - • • ,0„_i[n - 1], 6»„[n - 1] - a, 
6'„+i[n- I],-- - ,6»ArJn- 1]]'^ 

As in the previous step, the receiver compares the RSS 
function corresponding to 9"[n\ with the largest recorded RSS 
function and broadcasts a one-bit feedback signal Bp back 
to transmitters. Upon receiving ”1” from the receiver, the n- 
th transmitter updates its phase with the inversely adjusted 
value, while no action is taken for the cases of feedback being 
”0”. The feedback signal Bp is generated as in (j^ with 9'[n] 
replaced with 9"H and the procedures throughout Step 3 and 
Step 4 can be characterized by the following expression: 


Algorithm 1 Deterministic Bisection Search Algorithm 
(DBSA)_ 

Initialize parameters 
repeat 

k = 0 

for i ^ 1 to Ns do 

for 7 ^ 1 to I (360/a) — 11 do 
9'=9,[0]+a-j 
if RSS(6»'[0]) > RSS(0[O]) then 
6i[0] ^ 6i'[0] 
end if 
end for 

end for 

for z ^ 1 to Ns do 
k = k + 1 

9i [fc] 9i[k — 1] + a 
if RSS(6»[/c]) < RSS(6»[fc - 1]) then 
9i[k] ^ 9i[k - 1] - a 
if RSS(6'[A:]) > RSS(6'[A: - 1]) then 
9[k] ^ 9[k-l] 
end if 
end if 

end for 

a 

a ^ - 

untd stopping criteria reached 


half of its current value, and return to Step 3. For notational 
simplicity, the iteration index is reset to 1 and 0[O] is replaced 
with 9[Ns]. 

Stopping criterion: Given a threshold value for the RSS 
function, if the RSS function is greater than or equal to this 
target value at any step of the algorithm, the phase alignment 
process is said to be completed. The threshold value may be 
based on any statistical information and has to be available 
and known to the receiver. 

IV. Deterministic Quadratic Equation Search 
Algorithm 

Our proposed D-QESA algorithm undergoes Ng rounds of 
RSS function evaluations corresponding to the Ng distributed 
transmitters. During the z-th round, only the z-th transmitter 
transmits for 2 additional iterations. We will describe how 
these 2 additional RSS function evaluations can be used to 
improve phase alignment for different wireless environment 
settings in the following subsections. To begin with, note that 
the RSS function in (|^ can be rewritten as 


r 9'[n], Bp = l 

9[n] = l 0"[n], Bp = 0,Bp = l . (9) 

[ 9[n—l\, o.w. 

Once completed with the updating process, the algorithm 
returns to Step 3 and proceed to the (n-|- l)-th iteration where 
the phase of the (rz -|- l)-th transmitter can be adjusted, or, if 
n = N, the algorithm enters Step 5. 

Step 5- Parameter Adjustment: Upon entering Step 5, all 
transmitters scale down the phase adjustment parameter a to 


RSS(0i[rz], ■■■ ,9n^ [n]) = VP EiV a* 


= VP afcNe^'^^H + a,[n]e^-«-N 

k^i 

= VP\rt + ti\ 

= VP\I (|ri| -k \ti\cosl3if + {\ti\ sin/3Q^ 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


where n = and Pi is the 

phase angle between complex numbers and ti. Note that 
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the time indices of j3i, ri, and ti are omitted for notation sim¬ 
plicity. Furthermore, if one can only adjust the beamforming 
phase Oi[n] for each iteration, the optimal strategy would be 
aligning the direction of and ti, i.e., rotating 9i[n] to cancel 
out the phase angle Pi. This is the main idea behind D-QESA. 


A. D-QESA for static channels 

The pseudocode of the proposed algorithm can be found in 
Algorithm 2. For static channels, the channel coefficients do 
not change with time, i.e., hi = aiC^’^h Under this setting, we 
elaborate the proposed D-QESA in a step-by-step fashion as 
follows. 

Step 1- Initialization: The system is initialized 
by randomly generating the initial phase 0[O] = 

, 6 * AT Jo]} € In the beginning, we set 

9cur,i = ^*[ 0 ] where Ocur.i represents the current phase of the 
i-th distributed transmitter. 0temp represents the temporary 
value of phase. The parameter n is the time index with an 
initial value of 0. The phase adjustment parameters a and 
p are known to all transmitters. The angles a and p are 
initialized as a = tt and p = 7r/2. 

As mentioned earlier, the aim of D-QESA is to compute the 
value of Pi based on 2 additional RSS function evaluations and 
feed it back to the i-th transmitters to adjust its phase. In the 
following steps, we elaborate on how this is done based on 
only 2 additional RSS function evaluations in each round. 

Step 2- The phase rotation of a and rj: Denote the initial 
RSS function by Mi. We hrst rotate the phase of the i-th 
transmitter by a to obtain an updated RSS function value of 
M 2 . If M 2 > Ml, Bp is set to 1 and is fed back to distributed 
transmitters. The i-th transmitter then keeps this adjusted phase 
since the RSS function is improved. Otherwise, Bp = 0 is fed 
back and the i-th transmitter reverts back to its original phase. 
Note that this step not only improve the RSS function but also 
make Pi, to be solved later, falls in [— 7r/2,7r/2], the beneht 
of which will become clear later. 

Next, we obtain the second RSS function evaluations by 
rotating the phase of the i-th transmitter by p. Through this 
step, we obtain an updated RSS function value of M 3 at the 
receiver. 

Step 3- Calculate the value of Pi and |U|: Note that Mi 
represents the value of the initial RSS function value. M 2 
and M 3 were also obtained through 2 additional RSS function 
evaluations. Based on ( [T^ , we can obtain 

Ml = P -|-2|ri||U| cos (Pi)) 

Mi = P (|ri J -f -f 2|ri||U| cos (J* -f a)) (14) 

Mi = P (|ri J -f -f 2|ri||U| cos (-f p)) 


Algorithm 2 Deterministic Quadratic Equation Search Algo¬ 
rithm (D-QESA) 

Initialize parameters 

repeat 

for z ^ 1 to Ns do 
test •(— 0 
n •<— n -|- 1 


e-. 


temp ^ 9cur,i 


Ml i RSS£emp(^Zemp) 
RSS[n] ^ Ml 
n ^ n -I- 1 

9temp t 9cur.i -\-a 

M2 ^ RSSiemp(^Zemp) 
if M 2 > RSS[n - 1] then 
RSS[n] ^ M 2 


^cur,i ^ 

test ^ 1 


-f a 


receiver feeds back a 1 bit to transmitter 
else 

RSS[n] ^ RSS[n- 1] 

receiver feeds back a 0 bit to transmitter 

end if 

if n < 3N then 

n n -|- 1 

if test = 1 then 

9temp t 9cur,i ~i~ Tj O. 


else 

9: 


temp 


77 


end if 

Mu i RSSiemp(^Zemp) 

{Ml + Mi)/2 


Pi -tr- arctan 


{Mj - X 
.{Mi - x) 


\ti\ 


\ 


X — \ — 


{Mi — x)"^ 

(cos A) 2 


2P 


cur,r ^ ^cur.i 


else 


Pi •<— arccos 


Ml - Mi 


2\U\^/2{M^ + M|)P - 4\U\^P\ 

^cur,i ^cur,i 

end if 
end for 

until stopping criteria reached 


Then, 


Note that we describes 3 quadratic equations with 3 inde¬ 
pendent variable |ri|, \ti\ and Pi. Therefore, we can obtain 
\ri\, |U| and Pi by solving these equations. Specihcally, let 
X = {Ml M|)/2, we can obtain 


Mi - x = 2P|ri||U| cos {Pi -f tt) 


Mi - x = 2P|ri||U| cos (^Pi -f 0 

Mj-x _ + 2 ^ _ - sm{Pi) 

Mi — X cos {Pi -I- tt) — cos{Pi) 


= tan(/3i) 


( 16 ) 


x = P (|rip -P lUp) 


(15) 
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Now, we can obtain pi by 

It is clear that \ti\ = ai can be obtained by plugging ( [T^ 
into Specifically, we obtain 


H = \ 


X — 



(M| - xf 
{cos PiY 


2P 


(18) 


Step 4- With the information of |ti|, adjust the phase in 
the second round-rohin and after: During the first round- 
robin, we calculate the values of /3i and \ti\ for all i. Note 
that one can simply calculate the value of (3i and feed it back 
to adjust the phase of the i-th transmitter. Obtaining the value 
of \ti\ seems irrelevant for the current round. However, it is 
important to note that \ti\ equals the magnitude of the channel 
gain ai from the i-th distributed transmitter to the receiver 
and does not vary with time for a static channel. Furthermore, 
we run D-QESA in a round-robin fashion until the algorithm 
converges. By calculating and keeping the value of \ti\, we can 
reduce the 2 additional RSS function evaluations into only 1 
since now we only have two unknown independent variables 
and only Mi and M 2 are needed when we reaches the second 
round-robin and after. Specifically, from ([T^, we have 


Ml - M| = 4:P\ri\\ti\ cosiPi) 
Ml + Ml = 2P|r,p -f 2P|f,p 

Then, we can derive 


(19) 


Ml + Ml 

ip 

Plugging ( [T9] | into ( fTS] ), we obtain 

Ml - Ml = 4P 


- \t^\ 


Ml 


2P 


( 20 ) 


A/f2 \ 

^ |fi| cos(/3i) (21) 


Then, we can solve for /3i to obtain 

13, = arccos { - M| _ \ 

\2\U\^2{Ml + Ml)P - A\t,\^Pp 

This extra information of allows us to cut down the 
convergence time by about 1/3 after the second round-robin, 
which is a significant improvement. 

Once f3i is obtained, we can predict and adjust the phase at 
the Pth distributed transmitter accordingly. Ideally, one can 
feed back Pi directly to the Pth distributed transmitter to 
achieve perfect phase alignment from the viewpoint of the P 
th transmitter. However, if there is only limited bandwidth for 
the reverse feedback link Bp, proper quantization is necessary. 
For example, if 2 bits are available for the reverse feedback 
link. Pi can be quantized as {—Stt/S, —tt/S, tt/S, Stt/S}. With 
this feedback information, the Pth transmitter can decide 
whether to subtract either one of {—Stt/S, —tt/S, tt/S, Stt/S} 
to achieve a higher value of the RSS function. Note that 
if Pi is 0 or close to 0 , no information is fed back to 


distributed transmitters and the phase of the Pth transmitter 
is left unchanged. 

Stopping criterion: Given a threshold value for the RSS 
function, if the achieved value of the RSS function is greater 
than or equal to this target value at any step of the algorithm, 
the phase alignment process is said to be completed. The 
threshold value may be set based on the statistical information 
of channels. 


B. D-QESA-E for equal channel gains 

For the case where we have equal gains for all channels, 
i.e., ai = a, Vz, further modifications can be made to speed 
up our proposed D-QESA. We term the modified algorithm as 
Deterministic Quadratic Equation Search Algorithm - Equal 
channel gains (D-QESA-E). The pseudocode of the proposed 
algorithm can be found in Algorithm 3. Details of the algo¬ 
rithm is elaborated in the following steps. 

Step 1- Initialization: The system is initialized 
by randomly generating the initial phase 0[O] = 

{0i[O],-- - , 6 *AT Jo]} € In the beginning, we set 

(^cur,i = ^[ 0 ] where (lcur,i represent the current phase of the 
Pth transmitter. Again, the parameters a and p are initialized 
to a = TT and p = Tr/2. As in D-QESA, Pi is the value that 
we seek to calculate and feed back to the Pth transmitters 
for phase adjustment. The value of the RSS function for this 
initial phase is again denoted by Mi. 

Step 2- The first phase rotation of a and 77 : Here, the 
procedure is same as that in Algorithm 2. Eirst, we rotate the 
phase by a and use a one-bit feedback to determine whether 
the phase rotation is beneficial or not. If so, the updated phase 
is kept. Otherwise, the original phase is used. An updated RSS 
function value is obtained and denoted as M 2 . Then, we rotate 
the phase by p to obtain another RSS function evaluation M 3 . 
Note that this step is only done once. After the first round, 
only one additional RSS function evaluation is necessary for 
reasons that will become clear later. 

Step 3- Calculate the value of Pi and \ti\: The calculation 
of Pi and \ti\ is almost the same as before. There is, however, 
an important difference. Eor the case of equal channel gains, 
Oi = a, Vz. This implies that |p| = \aie^^'^‘^^ = a. After 
the first round of phase adjustment, we obtain |fi| = a and 
therefore, only two independent variables Pi and fp remain 
for the ensuing rounds. This important observation suggests 
that only Mi and M 2 are needed for the calculation of Pi. 
That is. 


3i = arccos 


Ml - Ml 


2\UW2{MI + Ml)P - ^U\^P^ 


(23) 


This modification allows us to cut down the convergence 
time of our proposed algorithm by almost 1 /3. This is a signifi¬ 
cant improvement for energy-efficient algorithms in distributed 
wireless sensor/relay networks. 

Stopping criterion: Given a threshold value for the RSS 
function, if the value of the RSS function is greater than or 
equal to this target value at any step of the algorithm, the 
phase alignment process is said to be completed. 
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Algorithm 3 Deterministic Quadratic Equation Search Algo¬ 
rithm - Equal channel gains (D-QESA-E) 


Initialize parameters 

repeat 

for i ^ 1 to Ns do 
test •(— 0 
n n + 1 

^temp ^ ^cur,i 

^^^tempi^temp') 

RSS[n] ^ Ml 
n ^ n + 1 

^temp ^ ^cur,i “t“ ^ 

M 2 RSSiemp(^iemp) 
if M 2 > RSS[n - 1] then 
RSS[n] ^ M 2 
^cur^i ^ ^cur,i + a 
test •(— 1 

receiver feeds back a 1 bit to transmitter 
else 

RSS[n] ^RSS[n-l] 

receiver feeds back a 0 bit to transmitter 

end if 

if n = 2 then 

n ^ n -f 1 

if test = 1 then 

^temp ^ ^cur,i “t“ ^ ^ 

else 

^temp ^cur,i “t“ ^ 

end if 


-^3 ^^^tempi^temp) 

X ^ {Ml + Ml)/2 


j3i •(— arctan 


\ti\ 


\ 


(M| — x) 
{Ml - x) 


X — \ x^ — 


{Ml - x)"^ 

(cos ft )2 


2P 


^ ^cur.i “f j 


-'cur,I ~ ^cur,i 


else 


f3i •(— arccos 


0. 


Ml - Ml 


^cur,z ^ 0cur,i “f 

end if 
end for 

nntil stopping criteria reached 


2|f,|V2(M2 + M|)P-4|f,|2p2 


C. Hybrid-QESA for time-varying channels 

In this subsection, we examine the RSS function of D-QESA 
under time-varying environments. Here, we assume that the 
time-varying channel phase is a one-step Markov process, i.e.. 


general. Therefore, a hybrid algorithm is proposed that com¬ 
bines the advantages of both kinds of adaptive beamforming 
algorithms. We term the proposed hybrid algorithm as hybrid- 
QESA. In hybrid-QESA, we initialize the algorithm by run¬ 
ning D-QESA for an entire round-robin, i.e., each distributed 
transmitter updates its own phase through D-QESA exactly 
once. Then, the algorithm switches to a random adaptive 
beamforming algorithm, BioRARSA2, which we proposed in 
our previous work |T5). Detailed definitions of parameters of 
BioRARSA2 can be found in Table 1 of p3] . The pseudocode 
of the resulting hybrid-QESA is described in Algorithm 4. 
The Hybrid-QESA combines the advantages of both rapid 
convergence of deterministic adaptive beamforming algorithms 
and the strong resistance to dramatic environments of random 
adaptive beamforming algorithms. This will be demonstrated 
in our numerical experiments. 


Algorithm 4 Hybrid-Quadratic Equation Search Algorithm 


Initialize parameters 

repeat 

for z ^ 1 to Ns do 
if n < “iN then 

test •(— 0 
n ^ n -f 1 

0temp ^ 0cur,i 

Ml i RSSiemp(^Zemp) 

RSSft] ^ Ml 
n ^ n 1 

0temp ^ 0cur,i ^ 

M 2 i RSSiemp(^Zemp) 
if M 2 > RSS[n - 1] then 
RSS[n] ^ M 2 

0cur,i ^ 0cur,i “f ^ 

test •(— 1 

receiver feeds back a 1 bit to transmitter 
else 

RSSft] ^ RSS[n- 1] 

receiver feeds back a 0 bit to transmitter 

end if 

n n 1 

if test = 1 then 

0temp ^ 0cur^i P 'h ^ 

else 

0temp ^ 0cur^i “f V 

end if 


M^3 i RSSiemp(^Zemp) 

x^ {Ml+Ml)/2 


ft ^ arctan 


0c 

else 




(M| — x) 

{Ml - x) 
+ Pi 


4>i[n] = (j)i[n - 1] ft [n] (24) 

where the sequence of ft[n] for i = 1,2,Ns consists of 
i.i.d. Gaussian noise. Namely, ft[n] ~ 7V(0,cr|). 

Eor time-varying channels, deterministic adaptive beam¬ 
forming algorithms are worse than their random counterpart in 


V. Performance analysis 

In this section, we analyze the convergence behavior of 
our proposed energy-efficient deterministic adaptive beam¬ 
forming algorithms, including DBSA, D-QESA, D-QESA-E, 
and hybrid-QESA. We are able to provide rigorous proofs of 
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for j ^ 1 to L.Helds do 

if RSSie^ur + ^) < RSS[n] then 
(5 =-(5 

end if 

if RSS(6'c„r + ^) > RSS[n] then 
Nj^ ^— 0 
uj„ <— 1 

repeat 

71 ^ n + 1 

^cur ^ ^cur “t“ ^ 

RSS[n] ^ RSS(0cuj') 

LOn ■<— LUn + 1 

untill RSS(0cur+^) < RSS[n] ora;„ > L.Swim 

else 

N'-p — Np -j-1 

end if 

i ^ i +1 

end for 

Avg. Swim ^ max(p, Y.]-L.Heid+i Tltids) 

k ^ k + l 

Afc ^ Afc_i- Avg. Swim 

if Np > Lp then 

A/j; i ^rst 

n ^ n + 1 

if RSS(0cMr) > RSS[77] then 
RSS[r7] ^ RSS[n- 1] 
else 

RSS[r7] •(— pt * RSS(0ciir) 

end if 

Np ■(— 0 

end if 
end if 
end for 

until stopping criteria reached 


convergence for all the above algorithms. Specifically, we can 
guarantee that all our proposed algorithms converge to the 
global optimal solution. 

A. Convergence analysis of DBS A 

The convergence behavior of DBSA is analyzed in the 
following Theorem. We are able to show that DBSA converges 
to the global optimal solution irrespective to the initialization 
parameters. 

Theorem 1: For the RSS function defined in 0, let 
{0 [t7]}^i be the sequence generated by D-QESA as described 
by Algorithm 1, where 6[n] = [0i[n], 92[n],--- , 6*ArJr7]]*. 
Then, the resulting sequence converges to the global optimal 
solution, i.e., lim„^oo RSS(6i[r7]) = RSS^aic = | Yh=i «i|- 

Proof: We design DBSA such that the phase for the z-th 
transmitter is only updated when the RSS function improves. 
This implies that the RSS function values achieved by 
the sequence are monotonically non-decreasing. 

Furthermore, it is clear that RSS(0[7i]) < RSSmax = | ^i\- 

Since RSS(0[77]) is upper bounded and monotonically non¬ 
decreasing, the convergence of DBSA is guaranteed 


by Monotone Convergence Theorem. However, this only 
guarantees that DBSA converges to local maxima. Fortunately, 
for the RSS function, all local maxima are global maxima 
(See 0 for details). Therefore, we can guarantee that DBSA 
converges to the global optimal solution. ■ 


B. Convergence analysis of D-QESA and D-QESA-E 


Here, we want to show that D-QESA and D-QESA-E indeed 
converge to the global optimal solution. Note that we do not 
need the property that all local maxima are global maxima as 
in the proof of DBSA. This makes our convergence analysis 
here applicable to a more general set of problem settings. 

Theorem 2: Eor the RSS function defined in 0, let 
{0[n]}^i be the sequence generated by D-QESA as described 
by Algorithm 2, where 6[n] = 92[n],"' j 

Then, the resulting sequence converges to the global optimal 
solution, i.e., lim„^oo RSS(6i[77]) = RSS^ax = | «i|- 

Proof: Erom ([T^, we can write the RSS function as 


RSS(0[n]) = 


Vp 


Ns 


VP\r,+f\ (25) 


By triangular inequality, we obtain 

RSS(0[n]) = VP\r,+f \ < /p(|r,| -f |f,|) (26) 


The aim of D-QESA is to align the phases of and f. 
Therefore, we obtain 


RSS(6>[n-f 1]) =/p(|r,|-f (27) 


This implies that RSS(0[n]) is a monotonically nondecreas¬ 
ing sequence with an upper bound of RSSmaa;- By Monotone 
Convergence Theorem, we can guarantee convergence of D- 
QESA. 

Next, we show that D-QESA indeed converges to the global 
maximum solution by contraction. Assume that D-QESA does 
not converge to the global maximum solution, there is at 
least one transmitter, say the z-th transmitter, not aligned with 
others. That is, 

RSS(6>[n]) = VP\r,+f \ < -f |f,|) (28) 

However, it is clear that RSS(0[n]) can be improved by 
running D-QESA for the z-th transmitter to align the phases 
of Ti and ti- This contradicts with our original assumption. 
Therefore, D-QESA will not stop at such a point, and it will 
stop only when 

RSS(0[n]) = Vp |r, + U\ = /p(|r,| + |f,|) (29) 

In other words, D-QESA is guaranteed to converge to the 
global optimal solution. ■ 

Eor the case of equal channel gains, we have the following 
theorem that describes and analyzes the convergence behavior 
of D-QESA-E. 

Theorem 3: Eor the RSS function defined in 
0 , let be the sequence generated by 







D-QESA-E as described by Algorithm 3, where 
6[n] — 6>ArJn]]*. Then, the resulting 

sequence converges to the global optimal solution, i.e., 
lim„^oo RSS(0[n]) = RSSmax = Nga. 

Proof: The proof is almost the same at that of Theorem 2 
and is not repeated here. ■ 

C. Convergence analysis of Hybrid-QESA 

The convergence analysis of Hybrid-QESA can be done 
through combining the convergence analysis of D-QESA and 
BioRASA2. The following theorem makes this statement more 
precise. 

Theorem 4: Eor the RSS function defined in (|^, let 
{0[n]}^i be the sequence generated by D-QESA as described 
by Algorithm 2, where 6[ri\ = 92[n],--- , 6*ArJn]]*. 

Then, the resulting sequence converges to the global optimal 
solution, i.e., lim„_j.oo RSS(6*[n]) = RSS,naa: = | Yh=i «i|- 

Proof: No matter how many rounds we run D-QESA before 
we switch into BioRASA2, the final phase configuration of D- 
QESA simply serves as a good initial point for BioRASA2. 
Since BioRASA2 is guaranteed to converge to the global 
optimal solution (See CD), Hybrid-QESA is also guaranteed 
to converges to the global optimal solution. ■ 

Remark 1: The proof of convergence for DBSA and 
Hybrid-QESA depends on the property that all local maxima 
are global maxima for the RSS function. However, we can 
guarantee convergence for D-QESA and D-QESA-E without 
this property. This indicates that convergence of D-QESA and 
D-QESA-E can be guaranteed for a wider set of problem 
settings, e.g., when the objective function is not the RSS 
function. 

VI. Simulation results 

In this section, we provide simulation results that evaluate 
the efficiency of our proposed deterministic adaptive beam¬ 
forming algorithms. We set the transmitted symbol power to 
be y/P — 1. In our simulations, all channels realizations are 
assumed to be zero-mean, unit variance i.i.d Rayleigh flat 
fading. All simulations are obtained with the same number 
of distributed transmitters, i.e., Ng = 100. Four types of 
channel and network settings are considered, i.e., noiseless 
channels, noisy channels, random node addition and removal, 
and noiseless time-varying channels. 

A. Noiseless channels 

In Figure 1 and Figure 2, we demonstrate how the the 
beamforming gain ratio evolves for noiseless channels. Each 
simulation curve is obtained by averaging over 100 randomly 
generated channel realizations. In Figure 1, we compare the 
convergence behavior between DBSA and D-QESA for dif¬ 
ferent number of feedback bits with the following schemes: a) 
one-bit random scheme proposed in Q, and b) BioRASA pro¬ 
posed in Q. We clearly observe that our proposed DBSA and 
D-QESA exhibit superior convergence behavior over random 
adaptive beamforming schemes. Comparison between DBSA 
and D-QESA reveals that better performance is achieved for 



Fig. 1. Comparison of the convergence behavior between D-QESA for 
different number of feedback bits (K = l, 25 inifite) with one-bit random 
scheme, BioRASA and DBSA for noiseless channels. 



Time slot 


Fig. 2. Evolution of D-QESA-E with different bandwidth of the reverse 
feedback link for the special case of equal channel gains, i.e., ai = a, Vi. 
Evolution of DBSA is also included for comparison. 


DBSA if there is only one-bit feedback. When more than 
2 bits are available for feedback, D-QESA becomes a more 
attractive alternative. We can observe that beamforming gain 
ratios for D-QESA almost all converge to a fixed value in 
about 300 iterations, which accounts for about one round- 
robin in D-QESA. That is, one round-robin is enough for the 
convergence of D-QESA. It is clear that faster convergence 
can be achieved with more number of feedback bits. This 
gain, however, diminishes quite fast as the number of feedback 
bits increases. For example, with 2 bits of feedback, 97% of 
RSSmaa: Can be reached in Figure 1. This indicates that D- 
QESA can also operate well with limited bandwidth of the 
reverse feedback link. 

In Figure 2, we compare the convergence behavior of D- 
QESA-E for different bandwidth of the reverse feedback link 
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Fig. 3. Comparison of the convergence behavior of D-QESA under different 
noise powers of —10 db and 0 db, respectively. The influence of D-QESA 
can be observed for noisy channels with one-bit and infinite number of bits 
of feedback . 


Fig. 4. Evolutions of D-QESA and DBSA are compare against the random 
one-bit scheme. The system with node addition probabilities ai‘e set to be 0.01 
and 0.05, respectively. 


with DBSA for = a. It is clear that beamforming gain ratios 
almost all converge to a fixed value in about 200 iterations, 
which accounts for about one round-robin in D-QESA-E. This 
demonstrates the rapid convergence of D-QESA-E. Indeed, 
there is about a 1/3 reduction in convergence time when 
compared with D-QESA, for which about 300 iterations are 
required before convergence. Eurthermore, D-QESA-E begins 
to outperform DBSA even with one-bit feedback before reach¬ 
ing its converged value of about 91% of RSSmaa;- 

B. Noisy channels 

In Eigure 3, we demonstrate the evolutions of D-QESA for 
noisy channels. Here, the performance curve is not obtained 
by averaging over 100 channel realizations since we want to 
observe the influence of noise. We compare the convergence 
behavior of D-QESA under different noise powers of —10 db 
and 0 db, respectively. When the noise power is same as the 
signal power, i.e., when noise power is 0 db, the beamforming 
gain ratio is poor and reaches only about 62% of RSSmax for 
one-bit feedback and 70% for infinite feedback bits. Eor a more 
practical noise level of —10 db, the beamforming gain ratio 
is a bit lower than that of a noiseless case and reaches about 
85% of RSSmax for one-bit feedback and 94% for infinite 
feedback bits. This demonstrates that D-QESA is still quite 
robust against noise for most practical scenarios. 

C. Random node addition/removal 

Here, we investigate the influence of random node addition 
and removal on adaptive beamforming algorithms. That is, we 
want to investigate whether our proposed algorithms are robust 
against uncertainties in wireless network topology. Note that 
we observe the RSS function value instead of the beamforming 
gain ratio here since random node addition and removal alter 
RSSmax such that the beamforming gain ratio may not be a 
good measure for performance. We choose two node addition 
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Fig. 5. Evolutions of D-QESA and DBSA are compare against the random 
one-bit scheme. The system with node elimination probabilities are set to be 
0.01 and 0.05, respectively. 


and removal probabilities of 0.01 and 0.05, respectively. In 
Eigure 4, we demonstrate the evolution of the RSS function 
under random node addition. We can see that fast initial 
convergence is maintained for both D-QESA and DBSA. As 
we move along further, the one-bit random scheme begins to 
outperform both D-QESA and DBSA, indicating that random 
adaptive beamforming algorithms are more robust against ran¬ 
dom node addition for wireless sensor/relay networks. Eor the 
case where the random node addition probability is set to be 
0.05, D-QESA outperforms DBSA slightly. This suggest that 
D-QESA might be more robust against network uncertainty. 

In Eigure 5, we demonstrate the evolution of the RSS 
function under different random node removal probabilities. 
A quite different phenomenon can be observed. It is clear 
that D-QESA and DBSA outperforms the one-bit scheme 
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Standard Deviation: 0.1 Standard Deviation: 0.15 



Fig. 6. Comparison of the beamforming gain ratio between Hybrid-QESA, 
D-QESA, DBSA, the one-bit random scheme, BioRASA, and BioRASA2 for 
different standard deviations (0.1 and 0.15) of the channel variation under 
time-varying channels. 


quite significantly. Furthermore, D-QESA demonstrates much 
superior performance than the other two algorithms, especially 
when the node removal probability is higher. Specifically, D- 
QESA is able to maintain its RSS function value to be above 
90% and 70% of RSSmaa: for removal probabilities of 0.01 
and 0.05, respectively. This is a surprising result since random 
adaptive beamforming algorithms are believed to be more 
robust against such uncertainties in general. This numerical 
experiment demonstrates that D-QESA is not only an attractive 
adaptive beamforming algorithm for static channels but also 
a robust algorithm against uncertainties in wireless network 
topology. This is a more desirable feature since random 
node removal is more detrimental for distributed beamforming 
schemes since there is a sudden decrease in the RSS function 
value. 


D. Time-varying channels 

In Eigure 6, we consider the case of time-varying channels 
described by (24 1 , where the standard deviation of is 
set to be 0.1 and 0.15, respectively. Here, we compare the 
convergence behavior of D-QESA, DBSA, and hybrid-QESA 
with the one-bit random scheme, BioRASA, and BioRASA2. 
When the standard deviation 5^. is set to be 0.1, we can see 
that the random adaptive beamforming algorithms have better 
performance than their deterministic counterparts. D-QESA 
still achieves a decent RSS function value but experiences 
severe fluctuation under time-varying channels. Hybrid-QESA 
combines the benefits of both types of adaptive beamforming 
algorithms and achieves the best performance over all the 
rest. Compared with BioRASA2, it can be observed that a 
fast initial convergence separates the performance of hybrid- 
QESA from BioRASA2. When % = 0.15, hybrid-QESA still 
outperforms most schemes. BioRASA2 has a slight edge in 
performance as we move further along. This suggests that 
random adaptive beamforming algorithms are still attractive 


alternatives under severe channel variations. However, strong 
performance for hybrid-QESA can be observed for both chan¬ 
nel settings. 


VII. Conclusion 

In this study, we proposed energy-efficient deterministic 
adaptive beamforming algorithms for distributed sensor/relay 
networks. Specifically, we proposed DBSA, D-QESA, D- 
QESA-E, and a hybrid algorithm, hybrid-QESA, that combines 
the benefits of both deterministic and random adaptive beam¬ 
forming algorithms. DBSA is inspired by the idea of bisection 
search. D-QESA is inspired by a shift of perspective that views 
each transmission as a RSS function evaluation and views 
the problem of adaptive beamforming for each distributed 
transmitter as solving quadratic equations with independent 
variables. We further provided rigorous convergence analysis 
for all our proposed algorithms and proved that the global 
optimal solution is reached for all our proposed algorithms. 
In our numerical experiments, we demonstrated that superior 
performance is achieved by our proposed DBSA and D-QESA 
over random adaptive beamforming algorithms for static 
channels. Surprisingly, D-QESA is also more robust against 
random node removal than random adaptive beamforming 
algorithms. Eor time-varying channels, hybrid-QESA indeed 
achieves the best performance since it combines the benefits of 
both types of adaptive beamforming algorithms. In summary, 
our proposed deterministic algorithms demonstrate superior 
performance both in terms of convergence time and robustness 
against channel and network uncertainties. This surprising new 
finding challenges the conventional belief that deterministic 
algorithms are faster in convergence while random algorithms 
are more robust against uncertainties. We hope that this work 
can generate more interests in the studies of both deterministic 
and random adaptive beamforming algorithms. 
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